arXiv:1503.07383vl [math.PR] 25 Mar 2015 


SINGULAR VALUES AND EVENNESS SYMMETRY IN RANDOM MATRIX 

THEORY 

FOLKMAR BORNEMANN AND PETER J. FORRESTER 


Abstract. Complex Hermitian random matrices with a unitary symmetry can be distinguished 
by a weight function. When this is even, it is a known result that the distribution of the singular 
values can be decomposed as the superposition of two independent eigenvalue sequences distributed 
according to particular matrix ensembles with chiral unitary symmetry. We give decompositions of 
the distribution of singular values, and the decimation of the singular values — whereby only even, 
or odd, labels are observed — for real symmetric random matrices with an orthogonal symmetry, 
and even weight. This requires further specifying the functional form of the weight to one of 
three types — Gauss, symmetric Jacobi or Cauchy. Inter-relations between gap probabilities with 
orthogonal and unitary symmetry follow as a corollary. The Gauss case has appeared in a recent 
work of Bornemann and La Croix. The Cauchy case, when appropriately specialised and upon 
stereographic projection, gives decompositions for the analogue of the singular values for the 
circular unitary and circular orthogonal ensembles. 


1. Introduction 


The ensembles of real symmetric random matrices OE„(r(;i) possessing an orthogonal symmetry, 
and complex Hermitian random matrices UE„(w 2 ) possessing a unitary symmetry, are specified by 
the eigenvalue densities 


(1) Pisixi,. ..jXn)= Cn^j3 ^^{Xk) ' |A(xi, . . .,Xn)\^ (/? = 1,2) 


with some normalization constant c„_/ 3 , each Xk restricted to the interval of support of wp{xk), and 
the Vandermonde determinant^ 


=det 






1 

6 


^ 2 ”' 


1 \ 

in 

C-7 


]J(6 - ij)- 

k>j 


Further, relating to a chiral unitary symmetry, there is the matrix ensemble chUE(r(; 2 ) with positive 
eigenvalues distributed according to the density, see [8, p. 717], 


(2) Pch{xi, • • • , Xn) — c„ 11 W2{xk) ' • • ■, . 


Key words and phrases, random matrices, evenness symmetry, singular values, gap probabilities. 
^Note that A(^i,... ,^n) ^ 0 if the arguments are increasingly ordered, ^ ^ ^n- 
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As in the theory of orthogonal polynomials, the w^^x) are referred to as weights. In fact the ensembles 
are often referred to by the name for the weights used in the theory of orthogonal polynomials. For 
example, OEAr(e““ is referred to as the Gaussian orthogonal ensemble. 

In this paper, as a unifying framework for examining eigenvalue properties under evenness 
symmetry, introduced into random matrix theory in the works [17, 8] and further explored in the 
Gaussian case in the recent works [7, 3], we study the structure of the singular values of ensembles 
OE„('u;i) with even weights supported on {—zu,zu) as given in Table 1. The ensemble of singular 

Table 1. admissible pairs of symmetric weights supported on {—zu,zu); a > —1 

case wi{x) W 2 (x) w 

Gauss e~^^ oo 

Jacobi (1 — x^Y (1 ~ 1 

Gauchy [l -p a;2^-(n+a+l)/2 _|_ 2,2^-(n+a) ^ 


values will be briefly denoted by |OE„(wi)|, in keeping with the relationship between the eigenvalues 
and singular values — since the ensembles are Hermitian, the singular values are the absolute value 
of the eigenvalues. Although defined according to the probability density function (1), we remark 
that each ensemble implied by Table 1 can be realised in terms of matrix ensembles defined by a 
distribution on the elements (see e.g. [10, Gh. 1-3]). 

Gentral to our discussion is the operation of decimation, which if applied to |OE„(wi)| results in 
the two ensembles 

even |OE„(wi)| and odd |OE„('u;i)|, 

where we define the even-location decimated ensemble even |OE„(r;i)| by taking the 2nd largest, 
4th largest etc. singular value, and similarly for odd |OE„(wi)|. The results will often depend on the 
parity /i of the underlying order n and we will, throughout this paper, write 

(3a) n = 2m + fj, (/i = 0,1), m = m + ^, 

that is, 

(3b) m = [n/2j , m = [n/2] , /x = [n/2] — [n/2j . 

Then, generalizing the corresponding result of Bornemann and La Groix [3, Thm. 1] for Gaussian 
ensembles, the following structure holds. 

Theorem 1. Let wp (jd = 1,2) be the weight pairs of the Gauss, symmetric Jacobi or Cauchy case 
as given in Table 1. Denoting equality of the joint distribution of two ensembles by =, there holds 

(4) even |OE„(t(;i)| = chUEm(a;^'^W 2 ) {n = 2m + fi). 

If we recall the superposition representation [8, Eq. (2.6)] 

(5) 1UE„(r; 2)| = chUEA(w2) U chUEm(a;^W2), 

of the singular values of the corresponding unitary ensemble UE„('u; 2 )) with both ensembles on 
the right drawn independently. Theorem 1 immediately implies the following remarkable relation 
between the singular values of OE(r;i) and UE(w 2 ): 
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Corollary 1. Let Wjs (f3 = 1,2) be the weight pairs of the Gauss, symmetric Jacobi or Cauchy case 
as given in Table 1. Then, with the ensembles on the right drawn independently, there holds 

(6) |UE„(w 2 )| = even|OE„(wi)| U even |OE„+i(wi)|. 

The superposition (6) bears a striking similarity with a corresponding superposition result for the 
eigenvalue distributions, see [12, pp. 185-186] or [10, §6.6], namely 

UE„(w 2 ) = even (OE„(rt;i) U OE„+i('u;i)). 

We proceed as follows: first, in Section 2 we give an overview of superposition and decimation 
results in random matrix theory known from previous studies, so as to properly set the scene for 
the present study and also as an opportunity to introduce the circular ensembles. In Section 3 
a factorized expression for the joint density of the singular values is obtained, with the proof of 
Theorem 1 in Section 5 following from this by integrating out the odd-location singular values. The 
success of this task is based on the notion of admissible symmetric weights, which we introduce in 
Section 4. There, Theorem 3 will give a complete classification of all the admissible weights, namely, 
they are exactly the Gauss, symmetric Jacobi and Cauchy weights (this is not to say that Theorem 1 
would not hold for other ensembles, but to point out that the method of proof is limited to those 
cases). In the first subsection of Section 6, some inter-relationships between gap probabilities are 
deduced from Theorem I. For n even, these have been obtained in the earlier study [8] without 
knowledge of Theorem I. We proceed to provide the necessary working to show that this is still 
possible for n odd. The relative complexity serves to further highlight the advantages of a viewpoint 
based on singular values. We conclude in Section 7 by presenting a number of new inter-relations 
between the spectra of circular ensembles, which follow upon the use a stereographic projection of 
the appropriate Cauchy weights to specify circular ensemble analogues of Theorem 1 and its various 
corollaries. 


2. INTER-RELATIONS KNOWN FROM PREVIOUS STUDIES 

2.1. Circular ensembles. Central to our theme is the operation of superposition, whereby eigen¬ 
value sequences from two independent ensembles with orthogonal symmetry are superimposed, and 
that of decimation, meaning in the present context that only those eigenvalues with a particular 
parity in the ordering are observed. The best known example of these operations involves not 
eigenvalues on the real line as in (1), but rather matrix ensembles with all eigenvalues on the unit 
circle in the complex plane. In fact such ensembles naturally follow from (1) with the Cauchy weight 

= (1 + ^ 2 )/ 3 (n-l)/ 2 -Ll • 

Thus, after making, for each eigenvalue, the change of variables 

(8) e*® = a; = tan(d/2), 

1 — ix 

corresponding to a stereographic mapping from the real line to the unit circle, one obtains the 
eigenvalue PDF on the unit circle 

(9) (x|A(e*®b...,e*^")|^ 

referred to, in the case /3 = 1, as the circular orthogonal ensemble COE„ and, in the case /3 = 2, as 
the circular unitary ensemble CUE„; see e.g. [10, Ch. 2]. 
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Let us superimpose two independent COE„ ensembles to obtain a new sequence of eigen-angles 

0 < 9l < 02 < ■ ■ ■ < 02n < 2Tr, 

and denote it by COE„UCOE„. It was conjectured by Dyson [6] and proved by Gunson [14] that 

(10) alt(COE„ U COE„) = CUE„, 

where the alt operation refers to the integration over alternate angles 01 , 03 ,, 02 n-i in the region 

& 2 j < (^ 2 j + l < ^ 2 j +2 (j = 0,..., n — 1) 

with 9q = 02 n — 274. 

The inter-relation (10) between eigenvalue distributions implies an inter-relation between condi¬ 
tioned gap probabilities. These are the probabilities, denoted by £’„_,g(A:; J; Wjs), or alternatively by 
En./sik; J;MEn,/ 3 iwi 3 )), that the matrix ensemble ME„_,g(w,g) contains exactly k eigenvalues in the 
interval J. Then as a direct combinatorial consequence of (10) one has [6, 15] (cf. also (15), (17)) 

(11) E„_2(fc;(-0,0);CUE„) 

n 

= {En,i{2{k - j); {-0, 0); COE„) + E^^,{2{k - j) - 1; {-0, 0); COE„)) 

j=o 

X (e„,i(2j; i-0, 0); COE„) + E„^,{2j + 1; (-0, 0); COE„)) . 

Closely related to the determinantal structure underlying the eigenvalue PDF (9) for the CUE„, 
together with the fact that this eigenvalue PDF is unchanged by complex conjugation, is the 
inter-relation [17] 

(12) |CUE„| = 0+(n-f 1) U 0"(n-f 1). 

As the name suggests, here 0^{n -I- 1) refers to the eigen-angles of matrices from the classical 
groups of the same name, chosen with Haar measure. Eigen-angles 0 and tt, which appear for purely 
algebraic reasons, are ignored and, since orthogonal matrices have real entries, for each eigen-angle 
0 7 ^ 0, TT, there is another eigen-angle —0, so that we take the one within the range 0 < 0 < tt only. 
The notation j • j now refers to the distribution of eigen-angles in the range 0 < 0 < tt, union the 
negative of the eigen-angles in the range —tt < 0 < 0. Though | • | has no effect on 0^(n -I- 1), this 
is not the case for the CUE„, where the eigenvalue distribution, and the distribution implied by 
|CUE„| are very different. 

As shown in [8, Eq. (2.6)], the analogue of (12) for Hermitian matrix ensembles with unitary 
symmetry is (5). In fact (12) can be deduced from (5) with the Cauchy weight W 2 {x) = (1 + a:^)“”, 
upon applying the change of variables (8) corresponding to a stereographic projection. On the RHS 
this requires the facts that under the change of variable x = tan6*/2 for each eigenvalue (see [8, 
Eqs. (2.24)-(2.28)], 

(13) chUEA((l+x2)-”) = Cl+(n-fl), chUE„(x^(l-h a:^)"") = 0-{n+l), 
on the LHS this change of variables simply gives 

|UE„((I+x2)-")| ^ |CUE„|. 
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2.2. Hermitian ensembles. Forrester and Rains [12] considered analogues of (10) for ensembles 
of Hermitian matrices. In keeping with above notations, let OE„(r(;i) U OE„(r(;i) denote the su¬ 
perimposing of two sequences of eigenvalues, independently drawn from OE„(u;i). Suppose the 
resulting eigenvalues are ordered Xi > X 2 > ■ ■ • > X 2 n, and let even (OE„(r<;i) U OE„(rt;i)) refer to 
the distribution of the even-location eigenvalues. We know from [12, pp. 186-187], see also [10, §6.6] 
that this is identically distributed to an ensemble with unitary symmetry, 

(14) even (OE„(wi) U OE„(wi)) = UE„(rt; 2 ), 

for the pairs (wi,W2) of weights given in Table 2 and furthermore, up to a linear fractional trans¬ 
formation, these pairs of weights are unique. The inter-relation between ensembles (14) has as an 

Table 2. pairs of weights satisfying Eq. (14); a > —1 

case wi(x) W2(x) support 

Laguerre (0,oo) 

Jacobi (1 — (1 — x)“ (0,1) 

immediate combinatorial consequence the inter-relation between gap probabilities 

2 k 

(15) En, 2 ik; (0, s); W 2 ) = ^ E„_i(2fc - j; (0, s); wi) (^E„,i(j; (0, s); wi) + En,i{j - 1; (0, s); wi)^ . 

3=0 

It is also fruitful to consider the superimposed and decimated ensemble even (OE„(/) U OE„+i(/)), 
thus involving one ensemble with n eigenvalues and the other with n -I- 1. It is shown in [12, pp. 185- 
186], see also [10, §6.6], that this, again, is identically distributed to an ensemble with unitary 
symmetry 

(16) even (OE„(wi) U OE„+i(wi)) = UE„(rt; 2 ), 

where (wi,W2) is any one of the pairs {wi,W2) of weights given in Table 3 (note that Table 1 gives 
the subset of even weights). As for Table 2 in relation to (14), these pairs of weights were shown to be 

Table 3. pairs of weights satisfying Eq. (16); a, 6 > — 1 


case 

wi(a;) 

W2{x) 

support 

Gauss 


e 

(— 00 , 00 ) 

Laguerre 

^(a-l)/2g-a;/2 

x^e-^ 

(0,oo) 

Jacobi 

(l + ^)(a-l)/2(l_3.)(&-l)/2 

{l+xY{l-xf 

(-1,1) 

Cauchy 

_|_ ^2^ —(n+a+l)/2 

(1 )-("+“) 

(— 00 , 00 ) 


unique up to linear transformation. An immediate combinatorial consequence for gap probabilities 
is the inter-relation 

2 k+l 

(17) En^2{k]Js\W2) = ^ En^i{2k + l- j-,Js;wi){En+i^i{j-,Js]Wi) +En+i^i{j - 1; J5;u;i)^, 

3=0 
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where Jg is a single interval either starting at the left boundary of support and finishing at s, or 
starting at s and finishing at the right boundary of support. 

Remark. Although it has no direct bearing on the present study, there is a decimation relation 
relating OE„(wi) for the weights in Table 3 to a corresponding PDF (1) with /3 = 4 [16, 12], which 
further generalises to a decimation relation reducing ensembles with /? = 2/(r+l),r€ TL'^ to 
ensembles with /3 = 2(r + 1) [9]. 


3. Joint Density of the Singular Values of Orthogonal Ensembles 


In this section we assume that Wi is an even weight function supported on the interval {—vj,vj). 
By symmetry, we can establish the joint density of the singular values by restricting ourselves to the 
cone of increasingly ordered singular values 

(18) 0 < CTi < • • • < CT„, 

this way parametrizing |OE„(uii)|. To simplify notation and to avoid case distinctions between odd 
and even order n in later parts of the paper, we introduce two further sets of coordinates for this 
cone. Writing, as detailed in (3), n = 2m + fj, and rh = m + n with /i = 0,1, the coordinates 

(19a) Xj = a 2 j-i (j = 1, ..., m), yj = a 2 j (j = 1,..., m) 

satisfy the interlacing property 

(19b) 0 < a;i < j/i ^ a;2 < j/2 < • • • ^ < 2/m < 

with formally adding, if /i = 1, the value ym+i = With x'^ and y^ denoting the x and y vectors with 
their components taken in the reverse order, so x'^ = {xm, Xm-i, ■ ■ ■ ,xi) and = {ym, y-m-i, ■ ■ ■ ,yi), 
we define, depending on the parity of n, the coordinates 

(20a) {t,s) = {y^,x^) (/a = 0), (t, s) = {x^,y^) iy = l), 


satisfying the interlacing property 

(20b) r<7 ^ ^ Si ^ ^2 ^ S2 ^ ^ ^ Sm ^ 0, 

again formally adding the value Sm+i = 0 if /i = 1. Since the mapping from cr = (cti, ... ,ct„) to 
either the pair of coordinates (x, y) or {t, s) is orthogonal, transforming the density between the 
three sets of coordinates is simply done by inserting new variable names for old ones. Note that the 
s variables parametrize the even-location decimated ensemble even |OE„('u;i)| while the /-variables 
do the same for odd |OE„(r(;i)|. We call them the even and odd singular values. 


By the evenness of wi the joint probability density of the singular values is, supported on (18), 


g((Ti,...,cr„) = n! ^ p(eiCTi,... ,e„CT„) = c„,in! • w((Tfe) • L>(cri,... ,CT„) 
ee{±l}" 


k=l 


with 




D(cti, ... ,(J„) = ^ |A(eicri,... ,e„ 

ee{±l}" 

Writing D(x; y) = D{ai ,..., cr„) in terms of {x, i/)-coordinates, Bornemann and La Croix [3, Eq. (11)] 
proved in two different ways the algebraic fact 

D{x\ 2 /) = 2” • A{xl, ...,x%)-yi--- ymA{yl ,..., 2 /^). 

Hence, we immediately get the following theorem: 
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Table 4. admissible symmetric weights wi(x) 


case 

parameter 

order 

wi{x) 

ZU 

Ok 

Pk 

p{x) 

Gauss 

— 

K < oo 

e-A/2 

OO 

1 

k- 1 

1 

Jacobi 

a > —1 

K < oo 

{1-x^y 

1 

1 

k-l 

1 — X' 

2cf< \ k 

2a -|- 1 -h A? 

Cauchy 

1 

“>-2 

K <2a 

(l-ba;2)-“"^ 

oo 

1 

k-l 

1 -1- a;' 

2a + 1 — fc 

2a -t 1 — fc 



^/W^(a + 1 ) 

2r(a+f) 
V^r(a+ I) 
2r(a+ 1) 


Theorem 2 . Letwi be an even weight on {—w,w). Then the joint probability density o/|OE„(wi)|, 
supported on the cone (19b), is given by 

( m \ / m 

]Jwi(a:fc) ■ A{xl,...,x%) I • I Y[ykWi{yk) ■ A{yl,... ,yl^) 

k=l / \fc=l 

with Cn = c„4n!2". 

Remark. Because of the interlacing in (19b), this factorization does not reveal an independence 
between the x and y variables. 


4. Admissible Symmetric Weights 

We call a smooth integrable weight wi : {—vo, zu) —>■ (0, oo) admissible of order n and mass 


( 22 ) 

if it satisfies the following properties 


29 = 




(a) wi is even; 

(b) Wi is normalized: Wi(0) = 1; 

(c) wi satisfies a three-term recurrence of antiderivatives of the form 

(23) J = -akx'"~^(l){x)wi{x) + Pk J ^'"~^wi{^)d^ {k = 1,2,..., k), 

with a smooth function p : {—w,w) —)■ (0, oo) and constants ak,Pk such that Pi = 0; 

(d) wi vanishes at the boundary: 

lim x^wi{x) = lim x^(j){x)wi{x) =0 (fc = 0,1,..., k). 

X—¥ZO X—¥-UJ 

Table 4 lists three cases of such admissible weights; by Theorem 3 below, these are all possible cases. 
By defining ao = 1, Po = 0 and 

1 r 

il){x) = --— -— / wi{(,)di {-w<x<w), 

(j){x)wi[x) Jo 
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the recurrence (c) extends to the case k 


0 if we replace a; ^ by By introducing the vectors 


<(x) 


X 




e M” 


y^u+2n-2 


) 


(u =-1,0,1), 


with the understanding that, instead of x the first entry of 7r’li(x) is in fact il){x), we can write 
the thus extended recurrence in compact matrix-vector form 


(24a) J wi(5)7r”(^) • wi(a;)7r”_i(a;) (u = 0,1, 2n + iy k + 2), 

(24b) Wi(x) = (l){x)wi{x), 


with a constant lower triangular matrix G having the numbers —a,y, —a^+ 2 , •.., —ai,+ 2 n -2 

along its main diagonal. In particular, there holds 

n—1 

(25) det = ( 1) Oi 2 k+iy 

k=0 

Since within the range of k restricted by the order k the constants ak and Pk given in Table 4 are 
strictly positive (with the exception of /3i = 0), we have An,i, > 0. We call wi = pwi the companion 
weight of wi{x) and observe that 

(26) lim wi{s)xl){s) = —9. 


In analogy to the results recalled in Section 2.2, we have the following uniqueness result. 

Theorem 3. Up to a rescaling of x, all possible admissible weights wi(x) are listed in Table 4- 
Actually, properties (b)-(d) of an admissible weight are sufficient for the conclusion to hold, that is, 
those properties already imply the evenness assumption (a). 


Proof. Let wi{x) be an admissible weight. Differentiating (23) yields 

(27) {x'^ - Pk + (k - l)ak4>)wi =-akxipwi)' (fc = 1, 2,..., k). 

Inserting a; = 0 gives 


Pk = {k- l)ak(j){0). 

Therefore, if Ofc = 0 for some positive integer k, we would get also that Pk = 0 and, hence, that 
rewiiOdf = o in contradiction to wi being positive. We conclude that 

akfO {k = 1,2 ,...,k). 

Inserting fc = 1 into (27) gives the differential equation 


(28) 


ipWiY = - XWi, 1111 ( 0 ) = 1. 

ai 


Inserting this expression for (pwi)' into (27) and rearranging, we get 


x^ - Pk . , x^ 


ak 


+ {k-l)P = 


ai 


{k = I,2,...,k). 
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Solving for (p gives 

1 Pk 
k-luk 

Since 4>{x) is assumed to be independent of k, we get that 
(29) 4){x) = 4){0) + Tx"^, ^ 0'2 ai 1 ttfc 

which can be solved for ak- 


1 Offc — «! 2 
h —f^ 

K — L akCXi 


Q;i 


Gfc = 


a2<Ti k — 1 akOti 
CXiOl2 


m- 

{k = 2,3,... ,k), 


a 2 + {k — 1)(q;i — q;2) 

Now, we distinguish four cases depending on whether ^(0) and r are zero or not. 

Case 1. (j){0) = 0 and r = 0, that is, (/> = 0. By (28) xwi = 0, which contradicts the positivity of wi. 

Case 2. (()(0) = 0 and r 7 ^ 0. By absorbing a rescaling of the ak into </> we can arrange for <j){x) = 

Now, solving the differential equation (28) for wi yields 

Wi{x) = CX 

with some constant c. For rui(O) = 1 to make sense, we would need the exponent to vanish, implying 
that already = 1. But such a weight would not satisfy wi(x) —)■ 0 as a; —>■ tu. 

Case 3. ^(0) ^ 0 and r = 0. By rescaling x we can arrange for ai = ±1 and (/> = 1. Now, solving the 
initial value problem (28) for rci yields 

w^{x) = e^^". 

From wi{x) ^ Q as x ^ w we get ai = 1 and w = 00 . This yields the Gauss case of Table 4. 

Case 4. 0(0) 7 ^ 0 and ai 7 ^ a 2 . By rescaling x and absorbing a rescaling of the ak into 0 we can 
arrange for 0(a:) = 1 ± Now, solving the initial value problem (28) for wi yields 

wi(a;) = (1 ± 

In the case (j){x) = 1 — we set a = — 1 + ^ and get, assuring integrability. 


wi(x) = (1 - x^Y 


a > — 1 , 


tZ7 = 1, 


1 


ai = 


2a + 2’ 

which yields the symmetric Jacobi case of Table 4. In the case 0(x) = 1 + we set a = 
get, once more assuring integrability. 


tY- and 

2a.\ 


Wi{x) = (I + X^) 


I 

a>- 2 , 


ai = 


1 

2a' 


which finally yields the Cauchy case of Table 4 (the only case where there is a restriction of the 
maximum order k that has to be checked). □ 


Remark. If 0(0) Y 0, Eqs. (28) and (29) imply that the logarithmic derivative of ihi = 4>wi, namely 

w[ X 

wi ai 0 (x)’ 

takes the form of a ratio of a linear and a quadratic polynomial. Hence, we immediately see that wi 
must be a classical weight. Because of the common denominator 0 in the logarithmic derivatives, 
the same conclusion holds for the weights iCi and W 2 = WiWi. Therefore, we could have finished the 
proof by checking properties (b)-(d) for each entry of a list of all classical weights. 
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5. Integrating Out the Odd and Even Singular Values 


5.1. Integrating out the odd singular values. We now prove Theorem 1. To begin with, we 
transform the joint density (21) to (s,t) coordinates, that is, 

(7(5, t) — Cn ' • ■ • j ' til—■ 5 

with functions 


(30) 






Likewise, we write for the same form of expression using the companion weight Wi instead of Wi. 

Now, Corollary 2 below shows that integrating out the odd singular values t subject to the 
interlacing (20b) gives the following marginal density of the even singular values: 

(31) Qeven('5l , . . . , Sjn ) — /j, ’ gfii^l-, ■ ■ ■ , Sm)gfi{^lj • ■ • j ^m) 

m 

= ■ p sl^W2{Sk) ■ A(s^, . . . , sjf (CC7 ^ Si ^ ^ Sm ^ 0) 

fc=l 

defining the associated weight function (cf. [12, Remark on p. 186]) 

(32) R'2 (s) = wi(s)wi(s) = ^(s)wi(s)^. 

Since the last expression in (31) is easily identified as the joint density of chUEm(a:^^W 2 ), see Eq. (2), 
we have finally proved Theorem 1. 

Remark. As a side product, the representation (31) shows that the normalization constant ^ of 
the joint density of OE(a;^^wi), if extended by symmetry to be supported on (0, oo)™, is given by 

ch _ A 

Lm,it — en,ir±rh.,l—fia . 

The integration is based on the following lemma and its first Corollary 2. 

Lemma 1. Let wi be the companion of the admissible weight wi. Then, there holds 

pX2 


d^i 


d^n det (u;i(^l)7r”(Cl) ••• 1Cl(Cn)’’‘"(?n)) 

-A det 


(u = 0,l). 


V 1 ••• 1 y 

Here, all integration bounds are within (O,^) and, in the case of a Cauchy weight, 2n + u ^ k + 2. 
Proof. Simplifying the notation to tTi,(x) = 7r”(x), we calculate by means of (24) and (25) 


nXn-\-l 


d^i--- I dfn det (wi(^i)7r^(^i) ••• wi(^„)7r^(5„)) 

•ai 2 

= det[ / wi(Ci)7ri/(6)c^Ci ••• 


P^n + l 


'^1 ('Cn)'^L' (^n)d<J?i 
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/ X2 


WlTT^-l 

• • WlTTiz-l 

V 

J 


— dct 


= (-l)"A„,„ 

( X2 
iTi(a;i)7r^-i(xi) witt^-i ••• witt^-i \ 

XI x„ 

1 0 ■■■ 0 J 

/wi(a;i)7r^_i(xi) ■u;i(a;2)7r^-i(x2) ••• ii;i(a;„+i)7r^_i(a;„+i)\ 

= An,i, det 

VI 1 ■■■ I ) 

In the last step we added the first column to the second, then the second to the third, etc. □ 
Corollary 2. Let g^, he as in (30) and put Sm = 0 if g = 1. Then, there holds 

(33) / dti I dt2 ’ ’ ’ / dt^ ? ^m) “ —/j, * ■ • ■ : (M “ 0, 1). 

9 Si J S2 9 

Here, all integration bounds are within (0,-ci7) and, in the case of a Cauchy weight, n = 2m + g ^ k + 2. 


(33) / dti / dt2" • 

J Si J S2 


Proof. Using the notation of Lemma 1, we hrst observe that 

(34) gf,{zi,. . . ,Zm) = det{wi{zm)T^Jf{z7n) ■■■ wi{zi)tt^{zi)) 

and the same for g^ with weight wi. Now, Lemma 1 yields, first using wi{s)tt™{s) —> 0 as s —>■ tu, 
that for /i = 0 


pSi 

dti / dt2 ■ 


— Am,i det 


dtm det{wi{tm)TrT{t-m) ■■■ Wl (fl)7r("(ti)) 

, /wi(Sm)7r™(Sm) ••• Wl(si)7r(("(si) 0 


1 IJ 

= ^m,idet(uii(s™)7r™(s™) ••• ?i;i(si)7r™(si)) 


and then, using 7r((‘(*‘^(0) = 0 and Wi{s)'iIj{s) —>■ — 0 as s — >■ tu, that for /x = 1 

pSi pSjn 

/ dti dt 2 --- dfm+i det (■u;i(fm+i)7r™+^(fm+i) • • • ■u;i(fi)7r™+^(fi)) 

J Si 9 S 2 9 Q 


— -^m+1.0 ■ d6t 


0 u;i(s„i)7r™j^^(s,„) ••• u;i(si)7r((“(''^(si) lim^^^ u;i(s)7r™('‘^(s) 


VI 1 ••• 1 1 ; 

= (- 1 )™A„+ 1,0 • det (u;i(sm)7r((‘j^^(sm) ••• wi(si)7r!(‘(''^(si) lim^^oo wi(s)7r™(''^(s)) 

, , J Wi{Sm)f’{s^) ••• Wi{si)->f{si) - 6 l\ 

= (-l) A„+i,o-det | ^ ^ ^ 


which finishes the proof. 


Wi(sm)7r™(sm) ••• i(;i(si)7r)"(si) 0 j 

= • det (wi(sm)7r("(sm) ••• wi(si)7r™(si)), 
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Remark. In the Jacobi case the multidimensional integral (33) can be recognized as a variant of the 
Dixon-Anderson integral [5, 1], well known in the theory of the Selberg integral, and also in the 
theory of /?-ensembles in random matrix theory [10, Section 4.2]. Specifically, in the statement of 
the Dixon-Anderson integral given in [10, Eq. (4.15)], cf. [5, Eq. ( 6 )], that is. 



j—l k—Q 


ntor(a.) 

r(Eto «0 


n 




valid for xq > Xi > • • • > x^h and aj >0 {j = 0,... ,m), we can reclaim the Jacobi case of (33) by 
the following substitutions of variables and choices of parameters: 

xo = 1, Xj = s], S,j=t) (j = 1,... ,to); 00 = 0 + 1 , Oj = 1 (j = l,...,m), 

and ara+i = 1/2, Sm+i = 0 if ^ = 1. 


5.2. Integrating out the even singular values. The following second corollary of Lemma 1 will 
allow us to integrate out the even singular values from the density (/(s; t). 

Corollary 3. Let g^, he as in (30) and put tm+i = 0 if g = 0. Then, there holds 


't2 


/■*™ , t 4 rl t ( 

ttSi • • • / dSm g^i[si, ...,Sm) = det 


dl — fj,{trh') 


di-tj.iti) , 


for g = 0,1 with 9 q(x) = 1 and 


di{x) = [ wi{^)d£,. 
Jo 


Here, all integration bounds are within (0,-ci7) and, in the case of a Cauchy weight, n = 2m + g ^ k + 2. 


Proof. Using (34) and Lemma 1 we obtain 
/ dsi * • * / ds^n g^{.S\, . . . , Sm) 

Jt2 Jt™+1 

= / dsi--- dsm det (wi(sm) 7 r™(sm) ••• R;i(si) 7 r™(si)) 

Jt2 Jim + l 

^l{tm+l)'!Tff_iitm+l) ••• Wi (t i )7r™_i (t i )^ 

which is already the assertion for g = 1. For g = 0, the assertion follows from further calculating 
^Wi{tm+i)n-i{tm+i) r(;i(tm)7r™i(im) ••• ^1 (U)7r!!*i (tl)\ 


det 


1 


1 


1 


0 Wi{tm)TrfJiitm) ••• Wi(ti)TT’fJ.^{ti) 

A 1 ••• 1 


= det 
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Table 5. companion weights wi(x) and integrals 0i(x) 


type 

wi(a:) 

wi(a;) 

6 * 1 ( 3 ;) 

Gauss 




Jacobi 


(1 - 

fl 3 ; 

Cauchy 

(l+a:2)-“-i 


/1 3 

a:-2Ti( 2’“+^’2’~ 


= (-l)'"det (wi(t^)7r™ (t„) ••• wi(ti)7r!!^^(ti)) 


= r-irdet( ••• 

[wi(t^)7r^ ^(im) ••• m(ti)7r^ ^(ti) 




-Wi(t 

m )^{t 

m) 




= det 

which finishes the proof. □ 

Now, by means of this corollary, the marginal density of the odd singular values is given as 

^odd (^17 • ■ • 7 ^m) — ' f?! —• ■ ■ j^l) ' det 



• U’l(6m)7r™^^ 

^1 — fii^tfi) 



Wl (tm)7r™ (t*) • • • Wi (tm)7r™ (tl 

7i-m(^a) ••• 7i-m(^i) 

(^m) * * * 


— Cn-^m,fi ' det 


• det 


with 


7^(x) = 0f,{x) = 


1 if /X = 0, 

[ if/x = l. 

Jo 


(n7 ^ ti ^ • • • ^ trfi ^ 0) 


Note that the two determinantal factors differ just in their last rows. It is this difference that prevents 
the expression from becoming a perfect square, which is in marked contrast with the marginal 
density (31) of the even singular values. 


6. Gap Probabilities 


6.1. A corollary of Theorem 1. Theorem 1 has an interesting implication in terms of gap 
probabilities, a notion that we recalled in Section 2. Specifically, we get the following result: 
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Theorem 4. The gap probabilities of the Gauss, symmetric Jacobi or Cauchy case of Table 1 of 
order n = 2m + /i satisfy 

En,i{2k + M - 1; {-s,s)]Wx) + En,i{2k + p; (-s, s); wi) = Em^ 2 {k-, { 0 ,s‘^)-,x^~^W 2 {x'^^‘^)Xx>o)- 


Proof. The change of variables Xk ^ Xk = y/xf, applied to the joint density Pch of the chiral ensemble 
chUE(x^''r(; 2 (a;)) yields 

Pch(a;i, ■ ■ • , Xm) dxi--- dXm = Pm, 2 {xi, ■■■, Xm) dxi ■ ■ ■ dXm 
where Pm ,2 is the density of \5E{x^~^W2{x^^‘^)Xx>o)- Hence, lifted to gap probabilities, we obtain 

(35) (0, s); chUE(x^^'u; 2 )) = Em^ik] (0, s^); x^^~^W 2 {x'^/'^)Xx>o) {h = 0, !)• 

By looking at pairs of consecutive values it is easy to see that the event that exactly k values of the 
decimated ensemble even |OE„(it;i)|, n = 2m + /i, are contained in (0, s) is given by the union of 
the events that exactly 2k + pL—1 or that exactly 2k + pL values of |OE„(wi)| are in that interval. 
Since these two events are mutually exclusive and since the singular values of OE„ contained in (0, s) 
correspond to the eigenvalues in (—s, s), we thus get from (4) and (35) proof of 

(36) En,i{2k + /I - 1; (-s, s);mi) + En,i{2k + (-s, s); wi) = Em, 2 {k; (0, s); chUE(x^^W 2 )) 

= Em,2{k] (0, s^); x^~^W2{x^^‘^)Xx>o), 


which finishes the proof. 


□ 


For even order (/i = 0), a first proof of this theorem was given by Forrester [8, Eq. (1.14)] using 
generating functions, Pfafhans and Fredholm determinants. For Gaussian weights, Bornemann and 
La Croix [3, Eq. (40)] recently settled the odd order case by using the more elementary techniques 
similar to this paper. 


6.2. Alternative derivation of Theorem 4. A natural question is to enquire if the proof of 
Theorem 4 for /i = 0 given in [8] can be extended to the case p, = 1. Here we will show that the 
answer is yes, although as is usual for methods based on Pfafhans in the study of random matrix 
ensembles with /? = 1 (see e.g. [10, Section 6.3.3]), the number of eigenvalues being odd adds to the 
complexity of the calculation. 

The hrst step is to introduce the generating function of the gap probabilities {En,p{k] J^wp)} 
according to 

OO 

En,0{J;£.;W0) = ^(1 - i)'^En,li{k-,J\Wp). 

The generating function can be expressed as the multidimensional integral (see e.g. [10, Prop. 8.1.2]) 

/ ZD pZD ^ 

dxi--- / dXn P[(l - f,Xxj(^j) -ppixi ,... ,X„). 

-ZD J — DD : -] 


( 37 ) 
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In terms of generating functions, the assertion of Theorem 4 in the case /i = 1 is equivalent to 

/ 1 ^2fc I Q 2 k +1 \ 

~ (2fc + l)!5C2'=+i j^2m+i,i((-s,s);e,m(x)) 

being valid for the weights in Table 1. It is this identity that we prove in the rest of the section. 

By making use of Pfaffians, (37) for /3 = 1 and R'i(a;) even can be expressed as a determinant. 

Lemma 2. Let Rj{x) he a polynomial of degree j for each j = 0,1,..., and furthermore require 
that Rj(x) he even (odd) for j even (odd). For wi{x) even we have 


(39) 

E 27 n+i,i{{-s,s);f]Wi) oc detv. 

where 


(40) 

V= ^[a2j-l,2fe] J = [&2j-l]j = l,...,m+l^ 

with 



(^3,k=7^ dxWi{x){l-fXxG{-s,s)) dy Wi{y){l - f,XyG{-s,s)) 
^ J —vj J —za 

X i?j_i(a;)sgn(y-x)i?fc-i(y), 

1 

(41) ^^^2] ^Xxe{-s,s))Rj-i{x)dx. 

The proportionality in (39) is such that the RHS is equal to unity when ^ = 0. 


Proof. Let h{x,y) = —h{y,x) and set 


^ f id{Xj,Xk)]j,k=l,...,2m+l [F{xj)]j=i^,,,^2m+l \ 

y {F{xk^]k—l,....2m+l 0 J 

It is well known (see e.g. [10, Eq. (6.81)]) that with h{x,y) = ^sgn.{y — x) and F{x) = ^ we have 

(43) PfX = 2-('"+i) Yl sgn{xk-x,), 

where Pf denotes the Pfaffian. Also, it is a simple corollary of the Vandermonde determinant identity 
that 

(44) det[i?fc-i(a;j)]j,fe=i,,„, 2 Tn+i oc (xk-Xj). 




Combining (43) and (44) shows 


\xk-Xj\ oc det[i?fc_i(a;j)]j_fe=i,„„2m+i Pf^- 


l^j<k^2m-\-l 
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The significance of the decomposition (45) for present purposes is that it implies a Pfaffian formula 
for the generating function E 2 m+i,i- Specifically, substituting the dehnition (1) of the joint density 
Pi{xi,... ,X2m+i) in (37) with /3 = 1, then substituting (45) we have 


dx 


2m+l 


/ Zu pZ 

dxi" ■ J 

2m+l 

^ (1 CXa:; e ( —s,s) )ll^l (^/) det X (Xj )] Pf ^ 

1=1 


oc Pf 


— [bk]k=l,...,2m+l 


0 


where aj^k, bj are given by (41), with the final line being a well known identity in random matrix 
theory [4], [10, Eq. (6.84)]. 

Finally, to obtain from this the determinant form (39), note that since (1 — ^Xxg{-s,s))wi{x) is 
even in x, and Rj{x) is even (odd) for j even (odd), we have that aj^k = 0 when j, k have the same 
parity, and bj = 0 for j even. Thus the nonzero entries in the Pfaffian (46) form a checkerboard 
pattern. Taking into consideration that aj^k is antisymmetric in the indices j, k, rearranging the 
rows reduces the RHS of (46) to 

^). 

V -Y^ 0^+1 J 

and this in turn is equal to det Y. □ 

At this stage the polynomials {Rj{x)}, apart from their degree and parity are arbitrary — a 
judicious choice takes us closer to establishing (38). For this, for a given W 2 {x) in Table 1, introduce 
the family of orthonormal polynomials {Pj(a;)}j=o,...,n such that 


(47) 


W 2 {x)pj{x)pk{x) dx = Sjk- 


In terms of these polynomials, and the pairs of weights as implied by Table 1, choose 
(48) i?o(a;) = l, R 2 j-i{x) =p 2 j-i{x), R 2 j{x) = - \—-^{^^^^^p 2 j-i{x) 


wi{x) dx \wi{x) 

where j = 1,2,.... The latter expression is even and a polynomial of degree j since 
1 d W2{x) X W2{x) 

wi{x) dx wi{x) «! ’ 


= (()(x) = even polynomial of degree 2, 


Wi{x)‘^ 

following from Eqs. (28) and (32), which are constitutive for Table 1. We then have, for j, fc = 1, 2 ... 


(49) 

and 


^1 = ^ y w^i(a;)(l - ^XxGi-s,s)) dx, b2j+i = 

^ pzz 

ai.2fe = X / (ixwi(x)(l - / dywiiy){l-^Xyei-s,s))sgii{y-x)p2k-iiy), 

^ J — zu J — zu 
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as well as 
(50) 02 ^ + 1 , 2 fe 


w ( 5 ) 

2?^7-TP2i-i(s) / Wi{y)p2k-iiy) dy 
Wl(S) Jg 

W 2 {y){i - ixy&{-s,s)fp 2 j-iiy)p 2 k-iiy) dy. 

-TU 


One immediate consequence of this choice is that it allows for a simple determination of the 
proportionality, 1/9 say, in (39). Thus with ^ = 0 we see that 62^+1 = 0 and 02 ^+ 1 , 2 fe = —dj,k, where 
to obtain the latter use has been made of (47). Consequently, cf. (22), 

(51) 6» = 6i|{=o = x / wi{x)dx. 

^ J — ZU 


Let C € 0{m) be a real orthogonal matrix, and define a set {< 72^-1 (a^)}i=i,...,m of polynomials by 



/ <li{x) \ 


/ Pi{x) \ 

(52) 

Qsix) 

= C 

Psix) 


\q2m-l{x)/ 


\P 2 m-lix)) 


If Y is defined as for Y but with each occurrence of p 2 j-i{x) replaced by q 2 j-i(x), we get 


(53) 


C 


Y 




= y. 


det Y = det F, 


where the latter follows from |detC| = 1. This allows us to make the same replacement in (48) and 
thus in (49) and (50) without effecting the representation (39) of the generating function. That this 
freedom leads to simplifications can be seen from the fact that {q 2 j-i{x)} remains an orthonormal 
set with respect to the inner product implied by (47), that is, 

/ ■cu 

W 2 {x)q 2 j-i{x)q 2 k-i{x) dx = 6jk, 

-ZU 

but can also be chosen to have an additional orthogonality as in the following lemma. 


Lemma 3. Define the projection kernel 

m 

(55) K{x,y) = {w 2 {x)w 2 {y)Y^‘^^P 2 k-i{x)p 2 k-i{y) 


together with the associated integral operator 

(56) Kf{x)=f K{x,y)f{y)dy (0 < s < tu). 

J — S 

This integral operator has eigenfunctions { 927-1 (a^)}i=i,...,m with the structure (52) for some real 
orthogonal matrix C, and furthermore 

(57) / W2{x)q2j-i{x)q2k-i{x)dx = i22j-i{s)6jk, 

J —S 

where 0 < ^ 27 - 1 ( 3 ) < 1 are the eigenvalues of K. 
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This functional analytic result is essentially due to Gaudin [13]; see also [10, p. 410]. The 
determinant of Y can be simplified by applying the elementary column operations of replacing 
column k for k = 1,..., n by column k minus 2 wi(x)q 2 k-i(x) dx times column n + 1. It is 
immediate that the entries in rows 2 ,..., n + 1 and columns 1 ,..., n are then given by 

/ ZU 

w 2 ( 2 /)(i - ^Xy&{-s,s)fq 2 j-iiy)q 2 k-iiy) dy 

-•07 

/ ZU pS 

W 2 (y)q 2 j-i(y)q 2 k-i{y)dy + {2s,-i'^) / W2(y)q2j-i(y)q2k-i{y)dy 

-zz j — s 

— djk (2^ ^ ^l'2j — ldjk — djk T (1 1) ^^2j — ldjkj 

where we have used (47) and (57) to obtain the last line. The entries in row 1, column 1, ..., n, after 
first simplifying the expression for ai^ 2 k in (49) by noting that the integral over y can be rewritten 
according to 


/ ZZ 

dywi(y){l - ^Xye(-s.s))sgn(j/ - x)p 2 k-iiy) 

-ZZ 

pZZ pZZ 

= ^Xx(^(-s,s) Wi(t)q2k-iit) dt + {1 - ^Xx(^{-s,s)) Wi{t)q 2 k-iit) dt, 

J S j X 


now read 


' dx 


\ dx 


pS pZZ 

(59) di. 2 fc=C(l-C)/ wi{t)q2k-i{t) dt / wi{x)dx + / wi(a;)(l - ^Xxe(-s.s))" 

J s J —S J —ZZ 

pZZ p'W p'W 

X / wi{t)q2k-i{t)dt - / wi{t)q2k-i{t) dt / Wi(a;)(l - ^Xxe(-s.s))' 

J X J S J —ZZ 

/ 'W pZZ pS pS 

wi{x)dx Wi{t)q 2 k-i{t) dt - {1 - - 1)'^) Wi{x) dx Wi{t)q 2 k-i{t) dt. 

-ZZ J X J —S J X 

The entries in the final column are unchanged by this process, and thus still have entries & 2 i-i ns 
specified in (49), with p2j-i{x) replaced by q2j-i(x). 

To summarize, we have shown with Eqs. (39), (51), (53), (58) and (59) that 


E27n+i,i{{-s,s);^]Wi) = -detF = 

C7 


1+^7 

_/+(!_ (^-1)2)^ eC3 


= det(/-(l-(e-l)^7?) + C 


cf + 4(l-(C-l)") 7 

_J+(1_(^_1)2)^ C3 


with D = diag(ui(s), 1 / 3 ( 5 ),..., U 2 m_i(s)), 7 a scalar and ci, 02,03 some column vectors with m 
entries that depend on s but not on The structure of the last formula is 

E 2 m+i,i{{-s, s)- roi) = E{1 - (^ - 1)^) + ^^(1 - (^ - 1)^), 


(60a) 
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where F{^) is a polynomial and 
(60b) 




Now, (60) is immediately amenable to the following simple lemma, which follows from direct 
compntation for the monomial basis {(1 — ^)'^}j=o,i, 2 ,...- 


Lemma 4. Let G'(^) be a polynomial. Then, for k = 0,1,2, 


/ 1 g2k ^ g2k+ 

~ (2fc +1)15^2'=+ 


^2fc+l 


^)G(l-(e-l)^) 


?=1 


(_l)fc gk 

k! 


kG(0 


«=1 


/ 1 g 2 k ^ 

~ (2k + l)!de^+^ 




J =1 


= 0 . 


Application of this lemma to (60) gives 

/ ^ g2k 2 ^2fc+l \ (~1)^ 

“ j2kTry.W^ jA 2 ™+i,i((-s,s);C;wi) = 

which finally proves (38) by the well known and readily established fact (see, e.g., [10, Exercises 9.6 
Q.3]) that 

E{0 = E^,2{iO,s^)-,^-,x^/^W2{x^/^)Xx>o)- 


7. Circular ensembles 

It was remarked in the paragraph including (13) that applying a stereographic projection to 
the eigenvalues in the appropriate Cauchy case of (5) gives (12). This transformation induces a 
natural definition of the decimated ensembles even |COE„| and odd |COE„|. Now, the analogue of 
Theorem 1 allows us to characterize not only the ensemble even |COE„| but also odd |COE„|: 

Theorem 5. Let fj, be defined as in (3) and, with sgn (x) = + for a; > 0 and sgn (x) = — for x < 0, 
define v = sgn (1/2 — pt). Then, the circular ensembles satisfy the inter-relations 

(61) even|COE„| = 0'^(n + l), 

(62) odd|COE„| = Q-'in + l), 

(63) |CUE„| = even|COE„| U odd|COE„|, 

where, in the last equation, both ensembles on the right are to be chosen independently. 

Remark. The last inter-relation should be contrasted with the trivial relation 

|COE„| = even|COE„| U odd|COE„| 

when both occurrences of COE„ on the right would represent one and the same ensemble instead of 
being independent. 
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Proof. The application of Theorem 1 to the Cauchy ensembles with weight (7) and a subsequent 
transformation to the circular ensembles by a stereographic projection of the eigenvalues transforms, 
by recalling (13), the inter-relation (4) into the hrst assertion (61). 

Next, we repeat these steps with the Cauchy weight 
(64) Wi{x)= (^^^2)(n-l+a)/2+l («>-!), 

which transforms by the stereographic projection into the circular Jacobi ensemble with parameter a 
[10, §3.9]. Though the resulting PDF becomes singular in the limit a —)■ — 1^, we know from working 
in the theory of the Selberg integral (see e.g. [10, Prop. 4.1.3]) that the limit effectively reduces the 
number of eigenvalues from n to n — 1 , by the mechanism of freezing one eigenvalue, taken to be 
at 0 = TT. This decouples but otherwise leaves the joint distribution of the remaining eigenvalues 
unchanged. Noting that the freezing of an eigenvalue at 0 = tt also has the consequence of replacing 
the even operation by the odd operation, and after applying analogous reasoning on the RHS of (4), 
we deduce the second assertion (62). 

Finally, by recalling (12), the last assertion (63) follows from (61) and (62). Alternatively, (63) could 
also have been deduced from ( 6 ) by the choice of the appropriate Cauchy weight, and appropriate 
interpretation of the weight in the second term as just discussed. □ 

Remark. Interestingly, the pathway to (62) via the limit a —)■ — I’*' in (4) with weight (64) can also 
be followed in the appropriate Laguerre and Jacobi cases of (16) to deduce (14). 

Analogous to the deduction of Theorem 4 from Theorem 1, as a corollary of Theorem 5, we get: 
Theorem 6. With fj, as in (3) and v = sgn(l/2 — /r), we have the gap probability inter-relations 

A„,i(2fc - 1 + m; {-0, Of COE„) + A„,i(2fc + /x; (-0, 0); COE„) = A„,2(fc; (0,0); 0+"(n + 1)), 

A„,i(2fc - (-0, e)- COE„) + E„,i(2fc + 1 - /x; (-0, 0); COE„) = (0,0); 0-"(n + 1)). 

In the case n even these inter-relations has previously been noted in [ 8 , Eq. (3.25)], where it is 
remarked that it allows the gap probabilities of COE„ to be expressed as simple linear combinations 
of the gap probabilities of 0^{n -\- I). One advantage of such expressions is that the ensembles 
0^(n -I- 1) are determinantal point processes (see e.g. [10, Ch. 5]), allowing the corresponding gap 
probabilities to be expressed as Eredholm determinants, which enjoy exponentially fast numerical 
approximation, and thus allowing for their efficient high precision computation [2]. Another advantage 
is that the gap probabilities for determinantal point processes can be shown to obey a local limit 
theorem in an appropriate asymptotic regime. The inter-relations then allow for the deduction of 
such asymptotic behaviour for the sum of neighbouring gap probabilities in COE„, for which no 
direct methods are known [ 11 ]. 

The gap probability inter-relation implied by (63) is exactly (11), even though the matrix ensemble 
inter-relation (10) used in its previous derivation is distinct from (63). This is a concrete example 
of the general fact that the family of gap-probability inter-relations specified by Theorem 4 or 
by Theorem 6 do not contain enough information to determine a particular matrix ensemble 
inter-relation, even though they are suggestive. 
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